On time evolution and causality of force-free black hole magnetospheres 



Amir Levinson 1,2 
ABSTRACT 

The rotational energy of a Kerr black hole is often invoked as the free energy source that 
powers compact astrophysical systems. An important question concerning the energy extraction 
mechanism is whether a Kerr black hole can communicate stresses to a distant load via a 
£\| . surrounding, force-free magnetosphere, and whether such structures are stable. In this paper we 

[ address this question. We first re-examine the properties of short wavelength force-free waves, 

D 1 and show that contrary to earlier claims, fast magnetosonic disturbances do affect the global 

electric current flowing in the system, even in flat spacetime; beyond the light cylinder the 
^sO ' perturbed charge density is on the order of the perturbed Goldreich-Julian charge density. This 

implies that the fast magnetosonic surface is a causal wind boundary. We then go on to study 
^vq , the evolution of a force-free system driven by a spinning-up black hole, by solving the time 

dependent Maxwell equations near the horizon using perturbation method. We find that the 
electromagnetic field, current and charge density grow exponentially in the linear regime. We 
conclude that after time t — 2Mln(l/2a 2 ), where a is the lapse function, the evolution will 
reach an adiabatic stage, whereby the solution is to a good approximation the static one, with 
the angular momentum of the hole being the adiabatic parameter. This implies that black hole 
magnetospheres are stable. A brief discussion on the relation between global current closure and 
' causality in force-free systems is given at the end. 
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Introduction 



Energy extraction from a Kerr black hole by a magnetized inflow is believed to be a plausible 
mechanism for powering compact astrophysical systems. The magnetic field lines threading the horizon 
may be anchored to a disk or a torus surrounding the black hole, where the extracted energy dissipates, 
or may extend to infinity, as envisioned in some applications of this mechanism to jet formation in AGNs, 
GRBs, and microquasars, in which case dissipation occurs at some distant load. 

Energy extraction along a magnetized flux tube that thread the horizon occurs when (i) the angular 
velocity associated with the flux tube is larger than zero and smaller than the hole angular velocity, and 
(ii) the Alfven point of the inflow is located inside the ergosphere (Takahashi et al. 1990). In the original 
model proposed by Blandford & Znjek (1977; hereafter BZ77), as well as in some other versions (e.g., van 
Putten 2001), it is conjectured that the hole rotational energy (or a fraction of it, as in van Putten's model) 
is extracted along an open, force-free flux tube that extends from the horizon to well beyond the outer 
light cylinder. In the context of ideal MHD it is clear that there must exist a region between the inner 
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and outer light cylinders where ideal MHD is violated, since the streamlines change direction (matter is 
inflowing into the horizon and outflowing to infinity). This region serves as a plasma source, e.g, via pair 
creation in a sparking gap. In the force free limit, in which inertia is neglected, this deviation is thought 
to be sufficiently small, to allow the force free approximation to be valid in the entire region between the 
horizon and the load - the region above the outer light cylinder, where the extracted energy dissipates (e.g., 
Blandford 2002). 

The question whether a force-free magnctosphere can exist and whether it is stable has been the subject 
of a recent debate (Punsly & Coroniti 1990, hereafter PC90; Beskin & Kuznetsova 2000; Komissarov 2001; 
Blandford 2001, 2002; Punsly 2003, van Putten & Levinson 2003). In particular it has been argued that a 
force free black hole magnctosphere is not a causal structure and, therefore, physically excluded. The point 
is that in the force free limit the electric current cannot flow across poloidal magnetic field lines, and is 
therefore conserved on magnetic flux surfaces. As a consequence the angular velocity, Q F , of a flux tube 
extending from the horizon is not a free parameter, but is determined by matching boundary conditions 
on the horizon and at infinity (see BZ77; Phinney 1983). This, according to PC90, violates the principle 
of MHD causality, because the inflowing magnetic wind must pass through the inner light cylinder before 
reaching the horizon and, therefore, cannot communicate with the plasma source region [e.g., the gap in the 
Blandford-Znajek model, or the torus in the model proposed by van Putten (2001; see also van Putten & 
Levinson 2003)]. PC90 concluded that the use of the Znajek frozen-in condition on the horizon to determine 
£If is unphysical, and that flF must be determined by the dissipative process that leads to ejection of 
plasma on magnetic field lines between the inner and outer Alfven points. 

Blandford (2001, 2002) claimed that a fast magnetosonic mode can propagate across magnetic field 
lines at the speed of light and carry information about the toroidal magnetic field and poloidal current to 
the plasma source even beyond the light cylinder. In the force free limit the fast critical surface approaches 
the horizon, and so stresses can be communicated to distant regions. Beskin & Kuznetsova (2000) reached 
similar conclusions by analyzing the Grad-Shafranov equation. A recent analysis of force free waves has 
been carried out by Punsly (2003), who has shown that linear perturbations of field aligned current and 
charge density cannot propagate along the fast wave characteristics. He therefore concluded that the use 
of the event horizon as a boundary surface in the force free limit is physically not allowed. The analysis of 
Punsly assumes the existence of a frame where the electric field vanishes and where Maxwell equations do 
not change form. We question this later assumption. In a frame rotating with the flux tube, for example, 
the electric field indeed vanishes. However, this frame is noninertial and Maxwell's equations should be 
properly transformed. In particular, the unperturbed charge density does not vanish in this frame and 
should be accounted for in the perturbed force-free condition for the waves. We discuss this point further in 
appendix A.l. Moreover, in Kerr spacetime proper account must be taken for the effect of frame dragging. 
Detailed investigation of force-free waves in Kerr spacetime is presented in Uchida (1997). Unfortunately, 
it does not address these issues in a clear way. Below, we re-analyze short wavelength disturbances and 
show that charge and current perturbations are in fact induced by fast magnetosonic disturbances, even 
in flat spacetime, and that beyond the light cylinder the perturbed charge density is of the order of the 
perturbed Goldreich-Julian (GJ) charge density. We also argue that frame dragging plays essentially a 
similar role as a unipolar inductor. This has already been pointed out earlier by Beskin & Kuznetsova 
(2000). To elucidate the effect of frame dragging on the time evolution of a black hole magnctosphere, we 
examine, in the second part of this paper, how a force free system evolves in the presence of a spinning-up 
black hole, by solving the time dependent Maxwell equations near the horizon using perturbation method. 
We find an exponential growth of toroidal magnetic field, current and charge density in the linear regime, 
and conclude that a black hole magnctosphere is stable. Our results appear to be consistent with recent 
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numerical simulations (Komissarov 2001; the rapid growth of the toroidal magnetic field found in the MHD 
simulations performed by Koide et al. [2002], and Koide [2003] is also consistent with the results described 
in sec. 2), and can also serve as a test case for such simulations. 



2. Analysis of force- free disturbances 

In the 3+1 formalism (Thorne, Price and Macdonald 1986; Beskin 1997) Maxwell's equations can 
be written, in terms of the electric and magnetic fields, and charge and current densities measured by a 
ZAMO, as 

1 <9B 

--+Vx(aE) = -(B.V/3)m, (1) 

1 (9E 1 

— + V x (aB) = -47raj + (E • V/?)m, (2) 

c at c 

V • B = 0, (3) 
V-E = 47rp e , (4) 

where the lapse function a, and the ZAMO angular velocity —/3, are given explicitly in appendix B. Here 
m = 94>4>^ < t > denotes the associated Killing vector field. In terms of the Boyer-Linduist coordinates it 
is given as m = Xsin6*/ p<f>, where the metric terms p and S are defined in appendix B. We denote the 
components of a vector in our local coordinate system by A = (A l7 A 2 , A 3 ), and choose A 3 to be the 
toroidal component (that is, in the direction of m). The force free condition takes the same form as in flat 
spacetime, viz., 

p e E+j x B/c = 0; E-B = 0. (5) 

From eq. (1) it is clear that since the ZAMO angular velocity is not conserved on magnetic flux surfaces, it 
introduces a coupling between the poloidal and toroidal components of the magnetic field. Furthermore, this 
term cannot be transformed away. The effect of this coupling on the evolution of the system is particularly 
important when the spacetime itself is changing with time. For instance, consider a situation wherein the 
angular momentum of the black hole undergoes low amplitude time variations, viz., a(t) = a + 5a(t), with 
5a « a. To first order in 5a, terms of the form — m(B • V)6fl(t), where B is the unperturbed magnetic 
field, and likewise for E, would appear in the wave equation for the perturbed electromagnetic field. These 
terms act as a driving force that excites waves, much like a Faraday disk with perturbed angular velocity, 
with one difference, that the driving force associated with frame dragging acts everywhere in space, not 
only on the horizon surface, although it declines steeply with radius. This qualitative discussion illustrates 
that spacetime does act like a unipolar inductor through the long range gravitomagnetic force, and that 
changes in spacetime are partly communicated to the plasma source by gravity. In the next section we 
present a quantitative treatment of this problem. 

The question whether the fast critical surface is a causal wind boundary still remains nonetheless, 
and is of relevancy also to super Alfvenic outflows in flat spacetime, as in the pulsar case. As we now 
show, if the unperturbed flux tube is rotating, then charge perturbations do propagate along cross-field 
fast characteristics even in flat spacetime. To elucidate this point we shall, in what follows, consider 
linear perturbations of a static solution of equations (l)-(5) far from the black hole. We can then set 
a = 1 and (3 = 0. Now, the unperturbed charge density approaches asymptotically the GJ value, that is, 
p e ~ pa j = E 2 cos9/2ttR, where R is a cylindrical radius measured with respect to the rotation axis and 9 
is the angle between the rotation axis and i\, and so in the force-free condition (5) the terms j x 5B and 



-4 - 



pSE are smaller by a factor of (kR)^ 1 than the terms Sj x B and Sp e E. This means that to zeroth order 
these terms can be neglected. As we shall show, the zeroth order solution gives the dispersion relations 
and eigen modes. However, it is clear that in order to maintain the asymptotic charge density near the GJ 
value, perturbations on the order of Sp e <~ 5E 2 /R are required; that is, the perturbed charge and current 
density appear only to first order. We shall therefore keep the terms associated with the unperturbed 
charge and current density in our analysis. We consider short-wavelength disturbances, and employ the 
WKB approximation, whereby the rapidly oscillating part of the perturbed quantities assumes the form 
expz[?/>(r) — ojt]. We define the wave vector of the disturbances as usual to be k = \7ip. We then obtain for 
the perturbed quantities: 

(fc 2 c 2 - uj 2 )SE - c 2 k(k • SE) = i4nuj6j, (6) 
ck x 5E = u;5B, (7) 
An6p e = i~k • SE, (8) 

cE5p e + cp e 5E + Sj x B + j x SB = 0, (9) 
B • SE + E ■ SB = 0. (10) 

These equations can be cast into the form Ay(k, u))5Ej = 0. The matrix Ay(k, ui) is derived for 
axisymmetric modes, viz., k$ = 0, in appendix A. In order to have a nontrivial solution, the determinant of 
Aij must vanish. The latter is given by eq. (A10). We find no corrections to first order to the dispersion 
relations of the intermediate (see eq [All)]), and fast (uj — kc) modes. In the remaining part of this section 
we analyze only the fast mode. The perturbed charge and current density of the fast mode vanish to zeroth 
order as inferred by Punsly (2003). However, using eq. (A14) we find that to first order the charge density 
is given by, 

k\\j p — kp e c 

where B p denotes the poloidal magnetic field, j p the poloidal current, and k\\ the component of the wave 
vector along B p . It is instructive to examine its behavior in the super- Alfvenic regime. Beyond the outer 
light cylinder the unperturbed poloidal current becomes purely convective, viz., j p = p e v p , and the toroidal 
magnetic field approaches B3 = E2 (e.g, Mestel 1999). Near the axis we therefore have B% ~ —VLpRBp/ c 
and k\\ ~ k\. Substituting these results into eq. (11) yields, 

^w^Sk^"^ (12) 

Thus, the perturbed charge density is approximately the perturbed GJ density beyond the light cylinder. 
This result follows essentially from the requirement on the wave to be force- free. Since u5p e = knSj p , it 
follows that perturbations near the fast critical surface will affect the global current flowing in the system 
and will be communicated to the source. The solution for the perturbed magnetic field can be found now 
from eq. (A13) and (7). In the region above the light cylinder it approaches 

SE 2 ~ -^■(l + k ]] /k)(n F R/c)5B 1 . (13) 

The analysis near the horizon is somewhat more involved, and will not be presented here. Recall, however, 
that the force-free condition (9) is valid also near the horizon in the ZAMO frame. By taking the projection 
of eq. (9) on the £3 direction we obtain: 



(p e -k-j)<y£3 + <5jxB = 0. 



(14) 



- 5 - 



Since SE3 7^ it is clear that current perturbations will be induced also in the vicinity of the horizon. 
However, in order to calculate the perturbed charge density, a solution for the eigen vector of the fast 
mode near the horizon must be found first. Our conclusion is that in cases in which the force-free limit is 
approached, the angular velocity of a rotating flux tube threading the event horizon of a Kerr black hole 
will be determined by the fast critical surfaces. 



3. Spacetime driven evolution of a force free magnetosphere 

Consider a black hole, initially non-rotating, embedded in an asymptotically uniform magnetic field 
directed along the rotation axis of the hole. We suppose that initially the magnetic field is described by the 
vacuum solution of eqs. (B6)-(B9), which is given by A t — F r g — 0, and 

= (B/2)r 2 sin 2 6>. (15) 

We further assume that the space surrounding the hole is filled with a tenuous plasma that maintains 
the system in a force-free state. Imagine now that the state of the black hole is perturbed adiabatically, 
such that it slowly acquires angular momentum, and that to a good approximation the spacetime can 
still be described by the Kerr metric with a = a(t). These temporal variations of spacetime will induce 
perturbations in the electromagnetic field. For simplicity we consider here the regime of slow rotation, viz., 
a/M « 1. For convenience we use the fields ip$, ip r , Bt, and defined in appendix B as our independent 
variables. To derive the equations for the perturbed electromagnetic field, we linearize eqs (B10)-(B17), 
using the hole angular momentum a as the small parameter. To zeroth order ipg = ip r = Bt = 0, and 
A<$, is given by eq. (15). The components of the metric tensor satisfy a 2 = 1 — 2M/r + 0(a 2 /M 2 ); 
A = r 2 - 2Mr + a 2 , p 2 = r 2 + a 2 cosd 2 ; and (3 = -2aM/r 3 + 0{a 3 /M 3 ). From eqs (B16) it is clear the is 
of order 0(a 2 /M 2 ), and that to second order in a/M the force-free conditions (B14) and (B15) reduce to, 

.,■ A .g g 

] =— — ] = -rcot0j , (16) 

-£±<f>,r 

Aip e = 4^Vv = r cot 6ip r - (17) 

From equations (16), (17), (BIO), (Bll) and (B13), we obtain a differential equation for the field ipg: 

(1 + A tan 2 0/r 2 )^ e , t ,t - -J—(Br,t),r + -^—^{B T ,t),e = 0, (18) 

r z sin r A cos 



with 



A sin 9 
Bta = 5 — 



/ a , x A , „ , , 6MBa . n 

(Awfl) r (tan c^ye) g 5 — sin cos 

r r z 



(19) 



From eq. (19) it is evident that the system is driven initially by the change of spacetime, specifically by 
the differential frame dragging term (i^F^g = —(6MBa/r 2 )sin8cos9. No evolution would ensue (to first 
order in a) in the absence of this term. In fact, it is easy to show that the coupling of the electromagnetic 
field to spacetime to first order in a is a consequence of the fact that the ZAMO rotation velocity — f3 is 
not conserved on magnetic flux surfaces. The above set of equations can be solved to yield the perturbed 
poloidal electric field and current, and toroidal magnetic fold, to second order in a/M. The charge density 
is then given by, 

^ + ^V),. + ^(sin0A*=O. (20) 
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Finally, linearizing eq (B12), we derive an inhomogeneous wave equation for the toroidal component of the 
perturbed vector potential SAa,: 



A 

SA 4>;tt + ~J 



sin 9 



,SA d 



= AnAj^, 



(21) 



with given by eq. (B16). At t = all perturbed quantities are zero, which defines the initial conditions. 

In the following we take for illustration a(r) = a [l — exp(— r/T)], where we define the dimensionless 
time r = t/2M. Near the horizon, the dimensionless variable x — a 2 /2 = A/2M 2 is small. Eq. (18) can be 
solved using perturbation approach. We suppose that the solution can be expanded as 



tpe = Y, n f n (r,x,6), 
where /„ is of order x n . To zeroth order we obtain 

d 2 f . Wa 



Or 2 /o AM 2 
B t .t = x2M 2 sin 9 (fo 



sm f cos ( 



35a \ 
AM 2 ) ■ 



The solution that satisfies the initial conditions fo — Bt — at r = is 



fo(r,x,9) = 



3Bao sin 29 



_ T 2 £ -(r/T) _ Tsinhr + cos hr + T 2 - 1 



8M 2 (T 2 - 1) 

There is no contribution to the electric currents j r and j to this order. To the next order we find 



fi(T,x,9) = f x 



3Bao sin 29 
8M 2 (AT 2 - 1) 



2T 2 (1 - e- T/I ) + -(cosh2r - 1) - Tsinh2r 

B T i = 2M 2 sin 9x [ (2/i - fox)dr' , 
Jo 



Ji 



3Ba cos 2 9 
AM 2 : 



(coshr — e T I T ) 



6T 



T 2 - 1 



T" 2 -I AT 2 -I 



(l-e- T/T ) 



(22) 

(23) 
(24) 

(25) 

(26) 
(27) 
(28) 



As seen, the perturbed force free field grows exponentially with an e-folding time 2M . It is clear that 
after time t ~ 2Mln(2M 2 /A) the term f\ becomes comparable to fo, and the above solution is no longer 
valid. We conjecture that after this time the evolution becomes adiabatic, in the sense that at any given 
time, r, the solution is given to a good approximation by the steady state solution with a = a(r) being 
the adiabatic parameter. Note that after time t ~ 2Mln(2M 2 /A) the toroidal magnetic field and charge 
density evolve close to their steady-state values, viz., Bt ~ (Qh / 2) F^g sin 6, and p e — a 2 j l ~ uB r cos 9, 
where — a/ 'AM 2 + 0(a 3 /M 3 ) is the angular velocity of the black hole and B r is the radial magnetic 
field. A more complete treatment will be presented elsewhere. The above suggests that the black hole 
magnetosphere is stable; the force-free plasma adjusts quickly to changes in spacetime. 



4. Discussion 



Based on the results found in sections 2 and 3, we argue that a force-free, black hole magnetosphere is 
globally causal, and that frame dragging plays an essential role in establishing such structures. We argue 
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that the force-free limit is a good approximation to ideal MHD systems when plasma inertia becomes 
negligible. Our analysis elucidates how the angular velocity of magnetic flux tubes penetrating the horizon 
responds to changes in spacetime, and confirms recent results found in numerical simulations (Komissarov 
2001; Koidc 2003). Moreover, the results of section 3 imply that a force-free magnetosphere is stable. 
Similar conclusions have been drawn by Beskin & Kuznetsova (2000), who present an elegant analysis 
of critical surfaces using the Grad-Shafranov equation. They show that the frozen-in condition on the 
horizon is an integral of the Grad-Shafranov equation, and contains no additional information, van Putten 
& Levinson confirm this point by showing that the frozen-in condition for a cold MHD inflow is merely 
a consequence of the fact that all trajectories approache those of a free-falling observer on the horizon. 
Whether the conditions required for a stationary force-free plasma to exist near a black hole are satisfied 
in nature is a different problem, that may involve issues concerning the microphysics of the plasma source, 
cross field diffusion, entrainment of matter by instabilities, etc. The other critical issue, yet to be resolved, 
concerns the global current closure in such systems. A particular example of a current circuit is discussed 
in van Putten & Levinson (2003). 

Even in flat spacetime one should be cautious in analyzing global causality. The analysis presented 
in Sec. 2 demonstrates that changes in the state of the system beyond the light cylinder (but not beyond 
the fast critical surface) can directly affect the global current flowing in the system. Moreover, current 
closure is a critical issue in such systems, and until it is properly modeled the question how the system 
adjusts to changes will remain unresolved. The wind region may be but one section of the global system; 
the return current may flow in subcritical regions, e.g., a cocoon surrounding a fast jet, and these regions 
may also affect the response of the system, regardless of our conclusion above. Recent analysis (Goodwin 
et al 2003) indeed suggests that boundary conditions beyond the light cylinder may considerably influence 
the properties of the wind in a pulsar. The response of the system may involve some feedback on the gap 
in relevant cases. In my view, sparking gaps as those envisioned in pulsar and black hole models are likely 
to have oscillatory behavior because steady state requires fine tuning of the microphysics. This may have 
important consequences for the properties of these systems that are yet to be explored. 

I thank J. Bekenstein, V. Beskin, Q., Luo, Y. Lyubarsky, D. Melrose and M. van Putten for useful 
discussions. This research was supported by an ISF grant for a Israli Center for High Energy Astrophysics. 

A. Solution of the wave equations 

Equations (7)- (9) can be solved for the perturbed current and magnetic field. We denote by B p and j p 
the unperturbed poloidal magnetic field and current density (Jz = j\ +if)> an d by 9j = j p /B p their ratio. 
Using the relations between the unperturbed quantities: E\ — (QfR/c)B 2 , and E 2 — — (ClpR/cjBi, and 
.73 = jpBs/Bp — p e E 2 /Bi, where fi_p is the angular velocity of the flux tube and R is a cylindrical radius, 
and the fact that the poloidal current flows along poloidal magnetic field lines, we can express the result 
in terms of the unperturbed magnetic field, charge density and poloidal current as: 5Bi = AijSEj, and 
Sji = DijSEj, where, 



Aij = (k 2 c/uj)S 3j , (Al) 
A 2j = -(k 1 /k 2 )A lj , (A2) 
A 3j = -{k 2 c/uj)5ij + (kic/u))5 2 j (A3) 
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_ Bioj 2 , , An A i4iruip e k 2 r , . 

^ = ~ k lBl + k2 B 2 k > + - kl B 1 + k 2 B 2 ^ 



u 2 



D 2j = -{k 1 /k 2 )D lj - —(kiSij + k 2 5 2j ), (A5) 

K 2 

D 3j = {B%/ B\){D\j - iA-KujOjAxj) + Mttw^Asj - <5 2jj o e /Bi) - (wfi F -R)(Mij + k 2 6 2j ). (A6) 

Substituting the above results into eq (6), and using eq. (10), yields a set of equations for the perturbed 
electric field: AijSEj = 0, where the matrix is given by 

Ay = Bj + (fifi2/w)(feiBi + k 2 B 2 )5 3j , (A7) 
A 2j = (fc 2 c 2 - w 2 )<5 23 - fc 2 fc 3 c 2 - D 2 j, (A8) 
Asj = -£>3, + (fc 2 c 2 - cj 2 )5 3j . (A9) 

Let us denote the poloidal magnetic field by B p , and the components of the wave vector parallel 
and perpendicular to the poloidal magnetic field by fc|| and k±, respectively. We can then write, 
k\B\ + fc 2 £? 2 = k\\B p , and k\B 2 — k 2 B\ = —k±B p . The determinant of Ay then takes the form, 



fcx(^ 2 -fc 2 c 2 ) 
k\\B p 



{(fc 2 c 2 - lu 2 )B 2 + (B 3 lu + n F Rk l{ B p ) 2 } . (A10) 



The requirement A = then yields the dispersion relations. We recover the familiar result; there are two 
modes, the intermediate mode satisfying, 

k\\cB n 



LU = 



Bl-Bl 



{-B 3 (n F R/c) ± [B 2 (l + n 2 F R 2 /c 2 ) - Bj} 1 ' 2 ) , (All) 



and having its group velocity directed along poloidal field lines, and the fast mode which satisfies ui = kc. 
There arc no higher order corrections to the dispersion relations. We now proceed with the analysis of the 
fast mode only. The solution for the eigen vector of the fast mode obtained upon substituting uo — kc into 
the equations for the perturbed electric field can be expressed now as, 

r F _ k2B 2 D 23 - k 2 (k 2 c 2 + D 22 )[B 3 + (np/^RkuBp} 

1 k 1 (k 2 c 2 + D 22 )[B 3 + (n F /u)Rk ll B p ]-k 2 B 1 D 23 2 ' ( ' 

6E = (k 1 B 2 -k 2 B 1 ){k 1 k 2 c 2 + D 21 ) 

3 fci(fcifc 2 c 2 + J D 21 )[S 3 + (n i r/a;) J RA;||S p ] - B X D 23 2 ' 1 ' 

Note that D 23 involves only first order terms. To zeroth order we therefore obtain 5E\ — —(k 2 /ki)SE 2 , 
implying that the perturbed charge and current density are indeed zero, as claimed by Punsly (2003). 
However, there are corrections to the next order. By employing eqs. (8), (A12) and (A13) one finds, 

Sp = k \pp - k P^ c 6E2 / A14 x 

kic[B 3 + (Q F /u>)Rk\\Bp\ - iAir(kiBi/k\\B p )c[k\\i P - kp e c] 

Note that the second term in the denominator is smaller by a factor of kR than the first term, and can be 
neglected to the lowest order. 



A.l. Wave equations in the rotating frame 



To clarify the discrepancy between our results and those of Punsly (2003), we note that Punsly assumes 
the existence of a frame where the electric field vanishes and where Maxwell's equations preserve their 
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form. While a local Lorentz frame where the electric field vanishes can be found (since B 2 — E 2 > 0), 
one should keep in mind that such a frame is only locally inertial, and there is no guarantee that the 
equations, when written in terms of the local quantities should preserve their form. To illustrate this 
point we now derive the wave equations in the rotating frame in the absence of a gravitational field. For 
convenience we use spherical coordinates, and denote quantities in the rotating frame by prime. In the 
non-rotating frame the line element is, ds 2 — —dt 2 + dr 2 + r 2 (d6 2 + sin 2 0d<j) 2 ). Performing a coordinate 
transformation to a frame rotating with angular velocity Q.f, viz., d<p' = d(j) + Qpdt (in units of c=l), one 
finds ds' 2 = -dt' 2 + dr' 2 + r' 2 d8' 2 + r' 2 sin 2 0'(d<j>' - tt F dt') 2 . All quantities can be transformed now to the 
rotating frame in the usual manner. It is readily seen that the only components of the electromagnetic 
tensor that are altered as a result of the transformation are those associated with the electric field: 
F[ 9 = F t e — VLfF^b, and F' tr = F tr — SlpF^. Since for a force free flux tube F t e — FIfF^b, and likewise for 
the radial component, we find that the electric field is indeed zero in the rotating frame. It is also easy to see 
that the only component of the 4-current that changes under the transformation is j'^ — + (r sin 8 ft _f)j', 
the other components remain unchanged. In particular the charge density does not transform away. Note 
that Gauss law is written in the rotating frame as, 

-1 [r 2 (F{ r + Q F f%.)] r - -^— [sm0(Fl r + Q F F^)] g = 4nj». (A15) 

Thus, the unperturbed charge density is given by the same equation as in the non-rotating frame, as 
expected from its transformation law. Amper's law becomes. 

{Ft + F F; r ), f + - JL^ ( S m6F^) e - 4nf r , (A16) 

and likewise for the 9 component. Linearizing these equations would yield the wave equations in the rotating 
frame where the electric field vanishes. Upon re-defining the fields Vv(e) = ^tr(e) + ^-^r(e)' ( w hich are 
of course just the components of the electric field in the non-rotating frame) one obtains exactly the same 
wave equations as in the non-rotating frame. 

B. Maxwell equations in Boyer-Lindquist coordinates 

The Kerr spacetime is a stationary, axisymmetric solution to the Einstein equations. In Boyer-Lindquist 
coordinates, (t, r, 6, <p), the line element can be written in the form, 

ds 2 = -a 2 dt 2 + g^dcj) + /3dt) 2 + g rr dr 2 + g ee d8 2 , (Bl) 

where 

p nr r, 2aMr p 2 2 o S 2 . 2n . 

a = -VA; = — ; g rr = — ; g ee = p ; g^ = u) = — sm 8 (B2) 

with A = r 2 + a 2 - 2Mr; p 2 = r 2 + a 2 cos 2 6; and S 2 = (r 2 + a 2 ) 2 - a 2 A sin 2 8. The parameter a = J/M 
represents the specific angular momentum. The units used above are geometrical units, in which r, M and 
a have units of length. 

The electromagnetic field in curved spacetime is describe by the tensor P" as usual. It obeys Maxwell's 
equations, 
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(B3) 
(B4) 



where j M is the generalized 4-current density. The force-free condition is F^j" = 0. In Kerr geometry 
Maxwell's equation reduce to the rather simple form: 



and 



Ftrfi + F r tt + F et ,r — F<ptM + Fg^ t — F^ t ,r + Fr<t>,t — 0, 



(B5) 



2 (^tr ~ flF^r) 



or 



■(Ftg — (3Fj,g) 



Asintf 



t<f> 



Asiu W-^) 



a" 



+ 



a z 94,4, 



' — ~2 {Ftr — 0F<t>r) 



a 
sin 9 

a 2 



+ 



(Fa - I3F^ 6 ) 



A sin 9 
Asinf 



-F 



r0 



-F 



r0 



(B6) 
(B7) 
(B8) 
(B9) 



It is convenient to define the quantities = —Eg/^A = T,(F t g — f3F,pg) / {p 2 A); 
ip r = -E r = T,(F tr - /3i^ r )/p 2 ; B T = (A sm9/p 2 )F r g; and j* = j* + /3f, where E r and Eg are 
the electric field measured by a ZAMO. The independent variables are then: tpg, ip r , Bt, and the toroidal 
component of the vector potential, A,/,, which in steady state defines magnetic flux surfaces. In terms of 
these quantities the inhomogeneous Maxwell equations reduce to the form 



(SW), t + ^ B T.fi = ^p 2 f, 



A sin 2 9 



A 



4>,t 



-Aa 



sm( 
1 /sinfl 



sinf 



The homogeneous Maxwell equation becomes 

„2 /f) 2 A 



Asin6> 



and the force free condition is written as: 



+ l3, r A<p,g - PfiAA,^ 



Af + A^gf + A^f - 0, 

A^rf + A^g f + A^.tf = 0, 



Vvj* - A^/f - 



A^gf - A^gf + 



p 



Asin6> 
A sin 9 ' 



B T f = 0, 



B T f = 0. 



(B10) 
(Bll) 
(B12) 

(B13) 

(B14) 
(B15) 

(B16) 
(B17) 
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